Condensation energy, charge ordering and inter-layer coupling in cuprates 
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The correlation between the condensation energy and the critical temperature is studied within a 
charge ordered superlattice bilayer model in which pairing is supported by inter-layer Coulomb en- 
ergy gain (potential energy driven superconductivity). The 2D pair-condensate can be characterized 
by a charge ordered state with a "checkerboard" like pattern seen by scanning tunneling microscopy. 
The drop of the c-axis dielectric screening can be the primary source of the condensation energy at 
optimal doping. We find that Coulomb energy gain occurs along the c-axis, which is proportional 
to the measured condensation energy (Uq) and to T c : E^ D » 2(f a f,/ao + i) 2 Uo ~ fcsTc and is due 
to inter-layer charge complementarity (charge asymmetry of the boson condensate) where £ a b is the 
coherence length of the condensate and ao m 3.9A is the in-plane lattice constant. The static c-axis 
dielectric constant e c is calculated for various cuprates and compared with the available experimen- 
tal data. PACS numbers: 74.20.-z, 74.25.-q, 74.72. -h 



Despite the enormous amount of effort spent in the last 
decades in resolving the fundamental physics behind high 
temperature superconductivity (HTSC) the mechanism 
is still puzzling [1]. It is more or less generally accepted 
now that the conventional electron-phonon pairing mech- 
anism cannot explain cuprate superconductivity, because 
as high a transition temperature as 1647^ (the record T c 
up to now [2]) cannot be explained by the energy scale 
of lattice vibrations without leading to lattice instability 
[1]. It is, however, already well established that much 
of the physics related to HTSC is in 2D nature, one of 
the basic questions to be answered still in the future is 
whether HTSC is a strictly 2D phenomenon or should be 
described by a 3D theory including all the 2D features of 
superconductivity. 

The systematic dependence of the transition tempera- 
ture T c on the c-axis structure and, in particular, on the 
number of CuC>2 planes in multilayer blocks are strongly 
in favour of the 3D character of HTSC. Variation of the 
thickness of ultrathin artificial HTSC compounds leads 
to the significant change of T c supporting the importance 
of IL coupling in cuprates [3] . It has also been shown re- 
cently that the c-axis transport of the hole-content to 
the doping site in the charge reservoir in the normal 
state (NS) of cuprates may be important in uderstanding 
HTSC [4]. Despite of the significant amount of published 
works which support the importance of the 3D character 
of HTSC 2D theories rule the literature in the last years 
[1]. This is partly due to the failure of the inter-layer 
(IL) tunneling model of Anderson [5,6] which was the 
prominant 3D theory of HTSC. Another reason could be 
that no sraigthforward dependence of T c can be seen on 
the IL distance and on other 3D parameters of cuprates. 
Therefore the 3D character of HTSC is more complicated 
and needs more theoretical and experimental work in the 
future. 



In this paper we propose a simple phenomenologi- 
cal model for explaining the 3D character of HTSC in 
cuprates supported by calculations. We would like to 
study the magnitude of direct Coulomb interaction be- 
tween charge ordered square superlattice layers as a pos- 
sible source of pairing interaction. Our intention is to 
understand HTSC within the context of an IL Coulomb- 
mediated mechanism. The IL charging energy we wish 
to calculate depends on the IL spacing (d), the IL di- 
electric constant e c , the hole content p and the size of 
the superlattice. A large body of experimental data are 
collected which support our model. A consistent picture 
is emerged on the basis of the careful analysis of this 
data set. Finally we calculate the static c-axis dielectric 
constant e c (that is normal to the a6-plane) for various 
cuprates which are compared with the available experi- 
mental observations. 

We propose to examine the following charge ordered 
superlattice model of pair condensation: A pair of charge 
carriers (2e) can in principle be distributed over 2/0.16 = 
12.5 CuOi unit cells in a square lattice layer if the 2e pair 
is composed of the optimal hole content po « O.16/Cm02 
at optimal doping [7]. However, allowing the phase sep- 
aration of hole-anti-hole pairs, every second unit cell is 
occupied by — 0.16e (anti-hole), and the rest is empty 
(holes, +0.16e), therefore we have 25 unit cells for a con- 
densed pair of charge carriers (FIG 1., that is the unit 
lattice of the pair condensate in the 5x5 supercell model) . 
Therefore, 2po = — 0.32e hole charge condenses to every 
second hole forming anti-hole sites. The 25 lattice sites 
provide the hole content of 25 x 0.1 6e = 4e and there- 
fore 2e excess charge in the anti-hole sites. The other 2e 
charge is used for neutralizing the holes at the anti-hole 
sites. In other words a Cooper wave-function can be dis- 
tributed on a 5 x 5 superlattice with a node of — 0.16e 
in every anti-hole site. The remarkable feature is that 




FIG. 1. The charge ordered state of the hole-anti-hole con- 
densate in the 4ao x 4ao charge ordered bilayer superlattice 
model. Note the charge asymmetry between the adjacent lay- 
ers. The bilayer can accomodate a pair of boson condensate 
(4e). The inter- layer charge complementarity of the charge or- 
dered state is crucial for getting inter-layer Coulomb energy 
gain. 



the size of the 5x5 condensate (four lattice spacings, 
4a « 15. 5A) is comparable with the measured small co- 
herence length £ a (, of single-layer cuprates (£ a & <~ 10 A to 
20 A) [8] . £ a b can directly be related to the characteristic 
size of the wave-pocket of the local Cooper pair (coher- 
ence area) [8]. The charge ordered superlattice bilayer 
(COSB) model can in principle be applied not only for 
po ~ 0.16e but also for the entire doping regime. 

IL Coulomb energy gain occurs only in that case when 
holes in one of the layers are in proximity with anti-holes 
in the other layer (FIG 1, IL electrostatic complemen- 
tarity, bilayer 5x5 (4a x 4a ) model). An important 
feature is then that the boson condensate can be de- 
scribed by an IL charge asymmetry. The IL coupling of 
the boson-boson pairs in the bilayer 5x5 model naturally 
suggests the effective mass of charge carriers m* « 4m e , 
as it was found by measurements [17]. The 5x5 model 
can be generalized to represent a N x N coherence area 
where N is the real space periodicity of the superlat- 
tice. In the N x N superlattice model the hole and 
anti-hole partial charges at each CuO-i lattice sites are 
q = ±4e/iV 2 [16]. Important to note that the charge 
sume rule holds for the characteristic bilayer with a co- 
herence area qf hole — 4e where q? hole is the partial 
anti-hole charge at the zth anti-hole lattice site. In other 
words a pair of a boson condensate can be accompanied 
within a COSB depicted in FIG 1. 

Important consequence of the model outlined in this 
paper is that the condensation energy (Uq) of the SC 
state can be calculated. We will show here that within 
our model the primary source of Uo is the net energy gain 
in IL Coulomb energy occurs due to the asymmetrical 
distribution of the condensed hole-charge in the adjacent 
layers (see FIG 1) below T c . 

Per definition Uo is the free energy difference of the 
normal state and the superconducting state, 



Tha a6-plane contribution to the SC condensation en- 
ergy AE ab is expected to be negligible at optimal dop- 
ing. That is because our Coulomb energy calculations 
indicate that the afe-plane Coulomb energy in the NS 
and in the SC state are nearly equal and therefore the 
ab Coulomb energy contribution to Uo A£° 6 « and 
also the kinetic energy of the pair-condensate Ej£ n is 
much smaller by several order of magnitude then IL cou- 
pling. The effective Coulomb interaction between spin- 
less point-charges may be approximated by the expres- 
sion V e ff(r) = e 2 /(4:TTe e c r), where e c takes into ac- 
count phenomenologically the dielectric screening effect 
of the IL dielectric medium and confined hole charge [16]. 
For instance we get for the prototypical 5x5 lattice 
~ -1.2 x 10~ 20 J afe-Coulomb energy both for the NS 
and for the SC COS using a simple point charge model 
for calculating the Coulomb energy. The details of this 
calculation can can found in ref. [16]. The calculated ab- 
kinetic energy of the 5x5 codensate is <~ 1.1 x 10~ 29 J 
if we use the very simple formula ("electron in a box") 
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where £ c is the c-axis coherence length. The IL Coulomb 
interaction (attraction) E* L w —1.3 x 10~ 21 J for the 
5x5 bilayer. It must be emphasized, however, that away 
from optimal doping the Coulomb energy of the conden- 
sate might affect the magnitude of the condensation en- 
ergy. The abrupt jump of the measured condensation 
energy Uo seen in the slightly overdoped regime [22] can 
be attributed to the increased AE ab contribution to Uo- 
In this article we restrict our study, however, to opti- 
mal doping and the study of the doping dependence of 
Uo is planned in the near future. The main contribution 
to Uq is then the IL energy gain AE^ L in the SC state 
can be given as follows, 
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where E^. L ' N and E^, L ' SC are the IL Coulomb energy in 
the NS and in SC state. We can further simplify Eq. (3) 
if the NS contribution to Eq. (3) is E* L ' NS « 0, which 
holds if IL coupling is screened effectively in the NS (large 
density of the hole content in the IL space, large e c (NS)). 
Again, if we assume e c « 33, we get E* L ' NS w +2 x 10~ 18 
J (12.3 eV) for a charge ordered 5x5 lattice which would 
be an extraordinarily large value for the IL Coulomb re- 
pulsion. Assuming, however, e c « 10000, we get the 
more realistic Coulomb repulsion of E[ L ' NS ~ +10~ 22 J. 
Indeed there are measurements for cuprates which indi- 
cate large e c in the NS in the range of 10 3 to 10 5 [27]. 
Extraordinarily large dielectric constant has also been 
found recently in perovskite materials [12]. Anyhow the 
"relaxation" of the huge IL Coulomb repulsion in hole- 
doped cuprates in the NS can not easily be understood 
without the consideration of a large e c . Nevertheless 
experimental measurements and theoretical speculations 



TABLE I. The calculated coherence length of the pair con- 
densate given in ao using the experimental condensation en- 
ergies of various cuprates and Eq. (6) at optimal doping. 
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ao Srf 3.88A, Uq is the measured condensation energy of various cuprates 
in fj, cV per unit cell at optimal doping. a from [28], b U ss 2 J/mol 
from [18 : 19], c from [8], d [30], e C/ ~ 12 - 16 mj/g from [20,31] and 
f„(, from [21], d (7 fs 11 J/mol from [22], c (7o S3 10 J/mol from [22], 
^ from [21], 9 from [22], % from recent measurements of Wang et al., 
tab ~ 23l(~ 5 - 6a ) [10], from STM images of ref. [14] £ ab RJ 4ao, h 
from [8], tab lc is calculated according to Eq. (6) and is also given in 
Table I and f;^ p is the measured in-plane coherence length given in 
ao ~ 3.9A. The notations are as follows for the compounds: Bi2201 is 
Bi 2 Sr 2 CuO e+s , LSCO (La 1 _ 85 Sro_i5Cu0 4 ), T12201 (Tl 2 Ba 2 CuO e ), 
Hgl201 (HgBa 2 CuO i+s ), YBCO (YBa 2 Cu0 7 ) and Bi2212 is 
Bi 2 Sr 2 CaCu 2 Og + s ■ 



suggest that e c in cuprates and in other materials with 
perovskite structure is strongly temperature and doping 
dependent [11]. The sharp decrease of e c below T c should 
also be explained within our bilayer model by the pair- 
condensation of the hole-content to the sheets. This is 
what can be seen in the c-axis optical spectra of cuprates 
[6,24]. The lack of the c-axis plasma edge in the NS is an 
obvious experimental evidence of the strong temperature 
dependence of e c . Then we have 
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where E^ L ' SC is the Coulomb energy gain in the SC state. 
Uq is the experimental condensation energy given per 
unit cell. Eq. (4) is generalized for multilayer cuprates in- 
troducing n. For single layer cuprates n — 0, for bilayers 
n = 1, etc. 

Inter-layer Coulomb coupling can be a natural source 
of the condensation energy and of HTSC. Checkerboard- 
like charge pattern (COS) seen experimentally [14,15] di- 
rectly leads to IL energy gain and to potential energy 
driven condensation energy if the hole-anti hole charge 
pattern is asymmetrycally condensed to the adjacent lay- 
ers (FIG 1). Assuming that thermical equilibrium occurs 
at T c for the competing charge ordered phases of the NS 
and the SC state the following equation for the conden- 
sation energy can be formulated using Eq. (4), 



2(n + l)N 2 U a « 2(n+l) 



u + 1 

a 



U nk B T c , (5) 



where kg is the Boltzmann constant. 

At a first look this formula seems to be unusual be- 
cause of the dependence of the condensation energy on 




2 4 6 

bilayer condensation energy 
FIG. 2. The critical temperature (K) at optimal doping as 
a function of the bilayer condensation energy (2(n+ l)N 2 Uo, 
meV) using Eq. (5). The straight line is a linear fit to the data. 
The slope of the linear fit is fcs which is a strong evidence of 
Eq. (5). 



T c . The available measurements of the condensation en- 
ergy on various cuprates show no correlation of Uq with 
the critical temperature. One of the important goals of 
this paper, however, to show that correlation can indeed 
be found with T c if T c is plotted against 2(n + 1)N 2 Uq 
(the bilayer condensation energy, FIG 2). In other words 
the condensation energy of the COSB system shows cor- 
relation with T c . In order to test the validity of Eq. (5) 
we estimate the coherence length of the pair condensate 
using Eqs. (5) and using only experimental data, 
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The results are given in Table I as £a& an d compared 
with the available measured £a& P - The agreement is excel- 
lent which strongly suggests that Eq. (6) should also work 
for other cuprates. The validity of Eq. (5) is clearly clar- 
ified in FIG 2 using mostly experimental information for 
C a6 , T c and U (values are given in Table I.). For T12201 
no measured £ a h is found in the literature, therefore the 
estimated £ a t, is used (given in Table I.). We will show 
that using this estimated value of N = 6(£ a {, w 5ao) we 
predict e c in nice agreement with optical measurements 
for T12201. Remarkable feature of FIG 2 is that the slope 
of the linear fit to the measured data points is precisely 
ks which supports Eq. (5). The underlying physics of 
HTSC seems to be reflected by Eq. (5): the equation cou- 
ples the basic observable quantities T c ,Uo and £ ab . We 
predict for the multilayer Hg-cuprates, HgT212 (T c = 126 
K) and Hgl223 (T c = 135 K) the condensation energies 
146 and 117 ^eV/u.c. w 4a [21]), respectively using 
Eq. (5). 

The expression Eq. (5) leads to the very simple formula 
for the critical temperature using Eq. (4) and a simple 
Coulomb expression for the IL coupling energy E I C L,SC 
[16] (T c » kz l El L > sc ) 



TABLE II. The calculated dielectric constant e c using 
Eq. (8) in various cuprates at the calculated coherence length 
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where £%% lc is the estimated in-plane coherence length given in ao ps 
3.9A unit, d is the CuOi plane to plane inter-layer distance in A, T c is 
the experime ntal critical temperature. e c is from Eq. (8). e^ xp are the 
measured values obtained from the following references: a [25], b [26], 
or from reflectivity measurements [30], cu p ps 55cm" 1 [33], A c ps 3^m 
[31], c from [30], d from reflectivity measurements: u p ps 60cm _1 [33], 
A c ps 0.9jum [13]. 
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is the inter-point charge distance and 

> diL, where diL is the IL distance (C11O2 plane 
to plane, i ^ j). n,m represent layer indexes (n ^ m). 
First the summation goes within the bilayer up to iV 2 
then the IL Coulomb interaction of the basal bilayer are 
calculated with other layers along the c-axis in both di- 
rection (n = 1,2) [16]. Ni is the number of layers along 
the c-axis. When Ni — > 00, bulk T c is calculated. T c can 
also be calculated for thin films when Ni is finite. e c can 
also be derived 



2 N, N 2 (n) (m) 
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where a c-axis average of e c is computed when Ni — > 00. 

The calculation of the c-axis dielectric constants e c 
might provide further evidences for Eq. (5) when com- 
pared with the measured values [11,26]. In Table II 
we have calculated the static dielectric function e c using 
Eq. (8) and compared with the experimental impedance 
measurements [27,26]. e c can also be extracted from 
the c-axis optical measurements using the relation [30] 
e c (w) = e c (oo)— c 2 / (u 2 A 2 ), where and uj p are the high- 
frequency dielectric constant and the plasma frequency, 
respectively [34]. c and A c are the speed of light and 

the c-axis penetration depth. At zero crossing e c (uS) = 
1/2 

and uj p — c/(A c e c (00)). Using this relation we predict 
for the single layer Hgl201 the low plasma frequency of 
uj p « 8 cm^ 1 using e c = e c (oo) w 6.5 (Table II) and 
A c « 8 /im [31]. 

If the physical picture derived from our model is cor- 
rect, it should be a guide for further experimental studies 
aiming to improve SC in cuprates or in other materials. 
This can be done by tuning the IL distance and e c in 
these materials. The application of this model to other 



class of HTSC materials, such as fullerides or MgB 2 is 
expected to be also effective. 
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